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Abstract 

This paper deals with the investigation of the computational solutions of an unified fractional reaction- 
diffusion equation, which is obtained from the standard diffusion equation by replacing the time derivative 
of first order by the generalized fractional time-derivative defined by Hilfer (2000), the space derivative of 
second order by the Riesz-Feller fractional derivative and adding the function (t>{x, t) which is a nonlinear 
function governing reaction. The solution is derived by the application of the Laplace and Fourier transforms 
in a compact and closed form in terms of the H-function. The main result obtained in this paper provides an 
elegant extension of the fundamental solution for the space-time fractional diffusion equation obtained earlier 
by Mainardi et al. (2001, 2005) and a result very recently given by Tomovski et al. (2011). Computational 
representation of the fundamental solution is also obtained explicitly. Fractional order moments of the 
distribution are deduced. At the end, mild extensions of the derived results associated with a finite number 
of Riesz-Feller space fractional derivatives are also discussed. 
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1. Introduction 



Standard (not fractional) reaction-diffusion equations are an important class of partial differential equa- 
tions to investigate nonlinear behavior of complex systems. Standard nonlinear reaction-diffusion equations 
can be simulated by numerical techniques, such as finite difference methods. The reaction-diffusion equation 
takes into account particle diffusion (diffusion constant and spatial Laplacian operator) and particle reaction 
(reaction constants and nonlinear reactive terms). Well-known special cases of such standard reaction- 
diffusion equations are the (i) Schloegl model, (ii) Fisher-Kolmogorov equation, (iii) real and complex 
Ginzburg-Landau equations, (iv) FitzHugh-Nagumo model, and (v) Gray-Scott model. These equations are 
known under their respective names both in the mathematical and physical literature. The nontrivial behav- 
ior of these equations arises from the competition between the reaction/relaxation and diffusion/transport. 
In recent years, interest is developed by several authors in the applications of reaction-diffusion models in 
pattern formation in physical and biological sciences and for describing non-Gaussian, non-Markovian, and 
non-Fickian phenomena in complex systems. In this connection, one can refer to Murray (2003), Kuramoto 
(2003), Wilhelmsson and Lazzaro (2001), and Hundsdorfer and Verwer (2003). These systems show that 
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diffusion can produce the spontaneous formation of spatio-temporal patterns. For details, see the work of 
Nicolis and Prigogine (1977) and Haken (2004). A general model for reaction-diffusion systems is investigated 
by Henry and Wearne (2000, 2002) and Henry et al. (2005), and Haubold et al. (2007, 2011, 2012). 

In this paper, we investigate the solution of an unified model of reaction-diffusion system (3.1) in which 

the two-parameter fractional derivative qDq'^ acts as a time-derivative and the Ricsz-FcUer derivative ^Dg 
as the space-derivative. This new model provides an extension of the models discussed earlier by Jesperson 
et al. (1999), Del-Castillo-Negrete et al. (2003), Mainardi et al. (2001, 2005), Kilbas et al. (2004), Haubold 
et al. (2007), Saxcna (2012), Saxena et al. (2010), Saxena et al. (2006 a,b,c,d), and Tomovski et al. (2001). 
The results are obtained in a compact form, which are suitable for numerical computation. Computational 
representation of the fundamental solution is explicitly derived and fractional order moments are investigated. 
For recent and related works on fractional kinetic equations and reaction-diffusion problems, one can refer 
to papers by Haubold and Mathai (1995, 2000) and Saxena et al. (2002, 2004a,b,c, 2006a,d, 2008). 

2. Unified fractional reaction-diffusion equations 

In this section, we will investigate the solution of the unified reaction-diffusion model (2.1). The main 
result is given in the form of the following 



Theorem 2.1. Consider an unified fractional reaction- diffusion model 

oD'^'-'Nix, t)=v xD^N{x, t) + cPix, t) 



(2.1) 



where r],t > 0,x G R;a,9,/3 are real parameters with the constraints < a < min(a, 2 — a), o-Df''' is the 
generalized Riemann-Liouville fractional derivative operator defined by (A7) with the initial conditions 



(l-i.)(l-^),0 



]N{x, 0+) = No{x); 0</i<l, Q<y<l 



(2.2) 



involving the Riemann-Liouville fractional integral of order (1 — — jj) evaluated for f — )• 0+ and 
\\ra\x\^rya N{x,t) = 0. Here xDg is the Riesz-Feller space- fractional derivative of order a and, asymmetry 
9 defined by (All), rj is a diffusion constant and (j>{x,t) is a nonlinear function belonging to the area of 
reaction- diffusion. Then for the solution of (2.1) subject to the above constraints, there holds the formula 



where 



a > and 



N{x,t)= f Gi{x-T,t)No{T)dT 
Jo 

/.t /.CO 

+ / (i-O""' / G2{x-r,t-0^{r,0dr 
Jo Jo 



Giix,t) = 



^H+v{l-IJ.)-l 



H. 



a\x\ 



2,1 
3,3 



(i,i),(;t.-KKi-M),^),(i,rt' 
(i,i),(i,i),(i,p) 



G2{x,t) 



a\x\ 



tH. 



2,1 
3,3 



(i,i),(M,^),(i,p)' 
(i,i),(i,i),(i,rt , 



(2.3) 
(2.4) 

(2.5) 



for a > where H^'\ is the familiar H-function [Mathai and Saxena (1978), Kilbas and Saigo (2004) md 
Mathai et al. (2010)], ^{jj) > 0, SR(;U + v{l - /x)) > 0. 

Proof: If we apply the Laplace transform (Erdelyi, et al. 1954) with respect to the time variable t and 
Fourier transform with respect to space variable x and use the initial conditions and the formulas (A8) and 
(All), then the given equation transforms into the form 



s''7V*(fc,s) -s"(''-i)7Vo*(fc) = -V i^i{k)N*{k,s) + (j>*{k,s), 



(2.6) 
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where according to the conventions followed, the symbol ~ will stand for the Laplace transform with respect 
to time variable t and * represents the Fourier transform with respect to space variable x. Solving for 
N*{k,s), it yields 



N*(k s) 

On taking the inverse Laplace transform of (2.7) and applying the formula 

o + 

where K(s) > 0, 3ft(a) > 0, SR(a - /3) > -1, it is seen that 

+ [ ci>*{k,t-o^^-'E^A-v€mndc 

Jo 

Taking the inverse Fourier transform of (2.9), we find 



(2.7) 



(2.^ 



(2.9) 



N{x,t) = 



2tt 



A^o*(fc)^^,^+.(i-M)(-»7 f'rocik)) eM-ikx)&k 



1 P 

^2^ Jo J * ~ eM-ikx)dk dC. (2.10) 

If we now apply the convolution theorem of the Fourier transform to (2.10) and make use of the following 
inverse Fourier transform formula (Haubold et al. 2007): 
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F-'[Ep,^{-rjt^i;i{ky,x] = ^<3^ 



(l,i),(7,^),(l,p)' 
(l,i),(l,l),(l,p) 



(2.11) 



where 5i(a) > 0, 5i(/3) > 0, 3?(7) > and p = {^^), it gives the desired solution in the form 

N{x, t)= [ Gi{x- T, t)No{T)dT 
Jo 



+ [\t-0^-' r G2{x-T,t-0<i>{r,0dTd^ 
Jo Jo 



where 



'0 Jo 



(2.12) 



Gi{x,t) = — J exp(-ifca;)i;^,^+^(i_^)(-?7t^V„(fc))dfc 



^fj.+v(i-ti)-i 



H: 



a\x\ 



2,1 
3,3 



(l,i),(,.+,.(l-Ai),f),(l,p)- 



(2.13) 



for a > and 



G2{x,t) = — I exp(-iA:a;)£;^,^(-77t^<(fc))dfc 



a\x 



tt2,1 
^3.3 



J. 



(i.^).(;',^),(i,p)' 



(2.14) 



where 3?(a) > 0, 3?(/u) > 0, 3?(/i + u{l — fj)) > 0. This completes the proof of the theorem. 
3. Special cases 

When = 0, the Riesz-Feller space derivative reduces to the Riesz space fractional derivative, and 
consequently we arrive at the following result recently given by Tomovski et al. (2011). 



Corollary 3.1. The solution of extended fractional reaction-diffusion equation 

oD^'^N{x,t) - r) ^D^N{x,t) = (j){x,t), x € R, t > 0, r? > 



with initial conditions 



< < 1, < M < 1 and < a < 2, lim N{x, t)=0 

x—>-zkoo 



]N{x,0+) = No{x) 



(3.1) 

(3.2) 
(3.3) 



where is a diffusion constant, oDq'" is the generalized Riemann-Liouville fractional derivative operator 

defined by (A7), xDq" is the Ricsz space fractional derivative operator defined by (A15) and 4>{x,t) is a 
nonlinear function belonging to the area of reaction-diffusion, is given by 



N{X, t) = r G^ix - T, t)No{T)dT 

Jo 



it 



Jo 



where 



for a > and 



G3{x,t) 



-H: 



a\x\ 



3,3 



Gi{x-T,t-^)<P{T,OdT d^, 

(1 i),(M+!^(i^p)4),(i,iy 
(li), (1,1), (1,1) 



G4(^,0 = ^^^3'i 



(i,i),(i,i),(i,i) 



, a > 0. 



(3.4) 
(3.5) 

(3.6) 



Note 3.1. Expressions for G3 and G4 can be obtained from Gi and G2 respectively by taking p = ^■ 

Similarly, if we set iVo(x) = (5(a;), </> = 0, where 5{x) is the Dirac delta function, then the theorem reduces 
to the following interesting result: 

Corollary 3.2. Consider the following reaction- diffusion model 



oD^'^N{x,t) = r] xD'^N{x,t), r] > 0, x e R, < a < 2, 



with the initial conditions 



.^(i-i/)(/*-i),Ox^. ^ 0+) = 5(x), < < 1, < < 1, lim N(x, t)=0 

+ X— >±oo 



(3.7) 
(3.8) 



where r] is a diffusion constant and 6{x) is the Dirac delta function. Then for the fundamental solution of 
(3.7) with the initial conditions (3.8), there holds the formula 



N{x,t) 



H. 



a\x\ 



2,1 
3,3 



(i,i),(^+..(i-/^),^),(i,p)- 
(i,i),(i,i),(i,rt 



(3.9) 



for a > where p = ^5-^ . 



Some interesting special cases of (3.9) arc enumerated below. 

(i) : When u = 1, then the operator D^'^ reduces to Caputo fractional derivative q defined in (A5) and 
the result (3.9) yields the fundamental solution of the space-time fractional diffusion equation obtained by 
Mainardi et al. (2001) and Mainardi et al. (2005). 

(ii) : Neutral fractional diffusion. We note that for p, = a, the solution (3.9) reduces to 



-H. 



a\x\ 



3,3 



{t r,^) 



(l,i),(c<+,/(l-a),l),(l,^))- 



(l,i),(l,l),(l,p) 



(3.10) 
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If wc further set i/ = I, then the time- fractional operator becomes the Capiito operator and neutral fractional 
diffusion occurs. It will be denoted by the conventional symbol A^^, (3.10) now simplifies into 



1 



7H, 



a\x\ 



1,1 
2,2 



" \x\ 









(3.11) 



which can be expressed in terms of a Mellin-Barnes type integral as 

/'7-t-zoo 



(3.12) 



{rj = 1, taken for simplicity). If the poles of the gamma functions occurring in the integrand of (3.12) are 
all simple, then evaluating the integral as a sum of the residues at the simple poles of T{s/a) at the points 
s = —an, n € No and r(l — s/a) at the points s = a + an, n € Nq we obtain the series representations 



and 



TTX — ' I 

00 

<(a;,t) = — ^sin[^(^-a)-""], 1< a; < 00. 



(3.13) 



(3.14) 



n=0 



Following the procedure adopted in Gorenflo and Mainardi (1997) and making use of the formula 

r exp(iQ;) 



^ r" sin(na) = exp(ma)] = ^{^^^^^X 



r sma 



for Irl < 1 



1 — 2r cos a + r"^ 

where a G R and it yields the interesting result given by Mainardi et al. (2001, 2005) 

1 y«-isin[f (a-fii)] 



K{x,t) 



iTT 1 + 2y° cos[|(q! - 6)] + y^a 



, ?/ = -,0 < a; < 00,0 < a < 2. 



(3.15) 



(3.16) 



Next, we derive some stable densities in terms of the H- functions as special cases of the solution (3.12). 

(iii): If we set /z = 1/ = 1, < a < 2, 6 < min{a, 2 — a} then (3.7) reduces to space-fractional diffusion 
equation, which we denote by L^{x), is the fundamental solution of the following space-time fractional 
diffusion model: 

^^^^ = ,D^N{x, t),rj>0,xeR, (3.17) 
ot 

with the initial conditions N{x,t = 0) = (5(a;), lim^^-too A^(a;, t) = 0, where r] is the diffusion constant and 
^(a;) is the Dirac delta function. Hence for the fundamental solution of (3.17) there holds the formula 



1 



L ^2,2 





(1,1), (i.p) ■ 


\x\ 


(i>i).(i,/'). 



, < a < 1, 1611 < a. 



(3.18) 



where p = ^2^- The density represented by the above expression is known as a-stable Levy density. By 
virtue of the H-function formula (Mathai et al., 2010) another form of this density is given by 



Liix) 



1 



Q!(r?t)< 



-H. 



1,1 
2,2 



" \x\ 


(o,i),(o,i)" 




(0,l),(0,p) 



, 1 < a < 2, \0\<2-a. 



(3.19) 



Remark 3.1. A general representation of all stable distributions in terms of special functions has been 
given by Schneider (1986). It is shown that the stable probability functions can be expressed by means of 
the H-functions, also see Uchaikin and Zolotarev (1999). We further note that Feller (1952) had derived 
the representations of the stable probability functions in terms of convergent and asymptotic power series 
in 1952. His results are revisited by Schneider, who had shown in case of L^{x) to restrict our attention to 
a; > 0, since the evaluation for a; < can be done by using the symmetry property L^{—x) — L^^{x). 

(iv): Next, if we take a = 2, 0</x<l, v = 1, = then we obtain the time-fractional diffusion, which is 
governed by the following time fractional diffusion model: 



N{x,t) = 7? —^N{x,t), 7? > 0, xeR, 



(3.20) 



with the initial conditions N{x, t = 0) = 6{x),\imx^±oo N{x, t) = where 77 is a diffusion constant and d{x) 
is the Dirac delta function, whose fundamental solution is given by the equation 



(v): Further, if we set a = 2,/i 
diffusion equation 



v = 1 and 9 



■^1,1 



" \x\ 






(1,1) . 



(3.21) 



d_ 
dt 



with initial conditions 



there holds the formula 
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N{x,t)=T] -^N{x,t), 
6{x) 



N{x,t = 0) 



lim N{x, t) = 0, 

X— >±oo 



then for the fundamental solution of the standard 

(3.22) 
(3.23) 



\x\ 






(1,1). 



(47r7?i)-5exp[-^, 



(3.24) 



which is the classical Gaussian density. For further details and importance of these special cases based on 
the Green function, one can refer to the papers by Mainardi et al. (2001, 2005). For f = the fractional 
derivative 0^0+" reduces to Riemann-Liouville fractional derivative operator q^D'^, defined by (A3), and the 
theorem yields 

Corollary 3.3. Consider an extended fractional reaction- diffusion model 

^''D'^N{x,t) = ri ,D^N{x,t) + cP{x,t) 



where ri,t > 0,x G R]a,9, fi are real parameters with the constraints 

< a < 2,|6'| < min(a,2 - a), 0<a<2 

where the Riemann-Liouville operator of order fi defined by (A3) has the initial conditions 

^^D^-'^N{x, 0) = No{x); ^^D^-^N{x, 0) = 0, < /x < 2, lim N{x, t) = 0. 

|x|^±C50 



(3.25) 
(3.26) 

(3.27) 



and is the Riesz-Feller space- fractional derivative or order a and asymmetry 9 defined by (All), rj is 
a diffusion constant and (l){x,t) is a nonlinear function. Then for the solution of (3.25) subject to the above 
constraints, there holds the formula 



N{x, t) 



+ 



G^{x-T, f)A^o(r)dr 

f{t - £,f-^ r Ge{x -T,t- O0(r, Odr d^ 
Jo Jo 



(3.28) 



where 



6 



a\x\ 



'3,3 



rjctc 



and 



1 rr2,l 

a\x\ '^'■^ 



(1 

(l,i),(A',^),(l,p)" 

(1 i),(i,i),(i,p) _ 



a > 



Ge{x,t) 

/ora > 0,3?(/x) >0,p+^. 

4. Fractional order moments 

In this section we will calculate the fractional order moments defined by 

/oo 
\xfN{x,t)dx 
-OO 

Using the result (3.9) and the fohowing definition of the Melhn transform 

POO 

M{f{t);s}= / t'-'f{t)dt 
Jo 



(3.29) 



(3.30) 



we find 



< x" > 



-f 

J — ( 



\x\^N{x,t)dx 



2iti+v('^-n)-i 



Jo 



(i,^),(/^+>^(i-;*),^),(i,p)' 



da:; 



The following formula gives the Mellin transform of the H-function (Mathai et al., 2010): 

+ BjS)}{UT=i r(l - aj - AjS)} 



f-OO 




(dp ,Ap ) 


/ x'-'H^f 

10 


ax 






(b,,B<,). 



Ax = a 



{nU+i r(i - - s,<5)}{m=„+i r(a,- + A^5)} 



(4.1) 

(4.2) 

(4.3) 
(4.4) 

(4.5) 



where -ro.va.i<j<rn^{^) < 5R(<5) < maxi<j<„ 3fi(^ji), |arga| < ^tt9, 9 = J^Jli^j - Z)^=„+i -Bj + 
J2^=i ~ J2'j=n+i -^3 ^ ^- ■'^PPlyiiig the above formula to evaluate the integral in (4.5), we see that 



< >= 2 A^„+.a-,H»(^)-i r(-A)r(i + j)r(i + A) 

T{-p5)V{i, + v{l-u) + lf)V{l + p6) 



(4.6) 



for -min{a,l} < sR((5) < 0. 

5. Computational representations of the solution (3.9) 

In this section we will derive the computational representation of the fundamental solution (3.9), which 
can be expressed in terms of the Mellin-Barnes type integral as 



N(x,t) = / J;/ ^ , — '—\ -^sm — (d-o)\\ Ms. 



(5.1) 



Let us assume that the poles of the gamma functions in the integrand of (5.1) are all simple. Now, evaluating 
the sum of residues in ascending powers of x by calculating the residues at the poles of r(l — s) at the points 
s = 1 + n,n € Nq and r(l — sa) at the points s = {1 + n) / a, n G Nq we obtain the following representation 
of the fundamental solution (3.12) in terms of two convergent series in ascending powers of x 



_ a;C.-i^^+.(i-M)-i ^ 

^„ r[M + - M) - (1 + n)M] 



r[l-a(l + n)] 



7r(77t'') 



(l+'«) ST,-X' 

sml^-: — -{9 - a)][ 



r;(l + /0/njr;l - (1 + »)/n] 



2a 

(1 + n)TT 

narj^t^ "'^t'" + - m) - (1 + n){ii/a)] 2a^ 

7 



■ sm 



r]o.to. 



where 



< 1. 



For v = 1 the time fractional derivative Dq^ becomes Caputo derivative qD^ and we obtain the 
following simplified form of the result given by Mainardi et al. (2001): 



^' TTivt^^) r[i - m(i ' ^ ^ '^^ ^ 



+ 



^ oo 



+ n)] 

r[(l + iij/n] 



2a 



(1 + Ti)7r — J 

7ra(r?^t^) n!r[l - (1 + n) Wa)] '^"^^^^^ " ""^^^^J 



(5.3) 



Further, if we calculate the residues at the poles of T{s) at the points s = —n, n £ Nq it gives 

N(x,t) = > — ^— 7sm — — (6'-a) 1< a; < oo. (5.4) 

^ ^ nx r(/i + z/(l - /i) + n/i) ^ 2a ^ '^^r]t^'^ ^ ' 

For u = 1 the above result simplifies to the following series representation: 

Ar/ ,N 1 ■^r(l + an) . ,^ .^^ a;" , „ ^ , ^, 

N(x,t) = — > -7- ^sm — — (6'-a) 1< a; < 00. (5.5) 

Finally, from (5.5) it follows that N{x,t) ^ jij for large \x\. 
6. Finite number of Riesz-Feller derivatives 

Following a similar procedure it is not difficult to establish the following 
Theorem 6.1. Consider an unified fractional reaction- diffusion model 

m 

oD^'-Nix, t) = Y^ Vj .D^l N{x, t) + (t>{x, t) (6.1) 

where i]j,t > 0,.x e R;aj,Oj,j = 1, m, /x, 1/ are real parameters with the constraints < aj < 2, \9j\ < 
min{Q;j,2 — aj}, oD^'''' is the generalized Riemann-Liouville fractional derivative operator defined by (A5) 
with the initial conditions 

j^a— )(i-A').Oj^(^^ 0+) = 7Vo(x); 0<Ai<l, 0<i^<l (6.2) 

involving the Riemann-Liouville fractional integral of order (1 — J^)(l — m) evaluated for i — >■ 0+ and 
\mY\x\^oo N{x,t) = 0. Here xDg^,,j = l,...,m are the Riesz-Feller space fractional derivatives of orders 
aj,j = l,...,m and asymmetry 6j,j = l,...,m respectively, defined by (A 10), 4>{x,t) is a nonlinear function. 
Then for the solution of (6.1), subject to the constraints above, there holds the formula 



4-/^+1^(1-^)-! 

N{x,t) = 



2-K 

r* 
10 



/OO ^ 
iV*(fc)i?^,^+.(i_^)(-t'^ ^r? {k)) eM-^kx)dk 
.7=1 



-1 pt poo ^ 

+ ^y^ e^-'y r{k,t-0E^A-t^J2vj i^%{k))eM-^kx)dk de (6.3) 



Some special cases of Theorem 6.1 are deduced below: 

(i): If we set 61 = ... = 6m = then by virtue of the identity (A13) we arrive at the following corollary 
associated with Riesz space fractional derivative: 



8 



Corollary 6.1. Consider the extended reaction- diffusion model 

m 

oDt'^Nix, t) = Y, Vj xDo' N{x, t) + 4>{x, t) (6.4) 

where rjj > 0, j = 1, m,t > 0, x G R, aj,j = 1, m, fi, v are real parameters with the constraints < aj < 
2, \9j\ < min(Q;j,2 — aj). qD^'" is the generalized Riemann-Liouville fractional derivative operator defined 
by (A7) with the initial conditions 

j^a— )(i-/^),Oj^^^^ 0+) = No{x); < /i < 1, 0<iy<l (6.5) 

involving the Riem,ann-Liouville fractional integral of order (1 — — /i) evaluated for t — > 0+ and 
lim|a;|_^oo A'"(a;, t) = 0. Then for the solution of (6.4) there holds the formula (6.3) with 'ipaj{k) replaced 
by \k\"i,j = l,...,m. 

(ii): If we further take = 1 in the above Corohary 6.1 then the operator [j^^^'^'^^^^)'*'] reduces to Caputo 
operator o^-Df defined in (A5) and we arrive at the following result: 

Corollary 6.2. Consider the extended reaction- diffusion model 

m 

^D!^N{x, t) = J2 .D^' N{x, t) + ^{x, t) (6.6) 

where all the quantities are as defined above, with the initial condition N{x,0-^-) = No[x); < aj < 
2, < n < 1, lim^x^^^ N{x,t) = 0. xDq' ,j = l,...,m are the Riesz space fractional derivatives of order 
aj,j = 1, ...,m defined by (AI4), ^{x,t) is a nonlinear function. Then for the solution of (6.6) there holds 
the formula 

1 /.oo m 

N{x,t) = — / N*{k)E^,r{-t''y^r)j \kfi)eM-ikx)dk 

+ ^y^C^-'y 0*(A;,t-Oi^/.,/.(-*''E^j |fcr)exp(-iMdfcdC (6.7) 



For m = 1 the result (6.7) reduces to one given by Tomovski ct al. (2011). 

(iii): If we set = then the Hilfcr (2000) fractional operator defined by (A7) reduces to Riemann-Liouville 
operator defined by (A3) and we arrive at the following 

Corollary 6.3. Consider an extended fractional reaction- diffusion model 

m 

^^D!^N{x,t) = .D;iN{x,t) + 4>{x,t) (6.8) 
i=i 

where the parameters and restrictions as defined before and with the initial conditions 

^^L)p^iV(x, 0) = Noix); I^^D'^~^N{x, 0) = 0, < < 2, lim N{x, t) = 0. (6.9) 

\x\—>oo 

Then for the solution of (6.8) there holds the formula: 

N{x,t) = -— / N*{k)E^^^{-t^^y^rji ^Pi^.{k))eM-^kx)dk 

J -00 

^2^/0 ^^'J •^*(^'*-^)'^''-''(-*''E^^- i^%)eM-^kx)dk de (6.10) 
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(iv): Finally, if wc set 9j ~ 0,j = 1, ...,m in Corollary 6.3 then the Riesz-Feller derivatives reduce to Riesz 

space fractional derivative and we arrive at the following 

Corollary 6.4. Consider an extended fractional reaction- diffusion model 

m 

^'^D>^N{x,t) = J2vj .D^'N{x,t) + cl>{x,t) (6.11) 

with the parameters and conditions on them as defined before and with the initial conditions as in ( 6. 7) then 
for the solution of (6.11) there holds the formula 

J./J-1 roo "1 

N{x,t) = -— / 7V*(fc)£^,^(-i^^77, |fcrOexp(-ifcar)dfc 

+ ^/^''"'/ <l^*{k,t-^)E^.^{-t<'y^rij\k\''^)eM-ikx)AkA^. (6.13) 

^TT Jo J-oo TTT 



7. Conclusions 

In this paper, the authors have presented an extension of the fundamental solution of space-time frac- 
tional diffusion given by Mainardi-Luchko-Pagnini (2001) by using the fractional order derivative operator 
defined by Hilfer (2000). The fundamental solution of the equation (2.1) is obtained in terms of H-function 
in closed and computable forms. Computational representations and fractional moments of the solutions are 
also obtained which will enhance the utility of the derived results in practical problems. Solutions of unified 
reaction-diffusion models associated with a finite number of Riesz-Feller space fractional derivatives are also 
investigated. 
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Appendix A. Mathematical preliminaries 
A generalization of the Mittag-LefHer function 

oo u 



^^(^) = E r(i + fca) ' ^gfi^W>o (^1) 

was introduced by Wiman (1905) in the generalized form 



EaA^) = J2 T{P + ak) ' ^^"^ > °' ^^^^ > ^^^^ 

The main results of these functions are available in the handbook of Erdelyi, et al. (1955, Section 18.1) and 
the monographs of Dzherbashyan (1966, 1993). The left-sided Riemann-Liouville fractional integral of order 
u is defined by Miller and Ross (1993, p.45), Samko et al. (1990), Kilbas et al. (2006) as 

1 /■* 

^^D^''N{x, t) = I^N{x, t) = —— {t- uY-^N{x, u)Au, t > 0, 3?(i/) > 0. (A3) 

r(z^) ^0 

The left-sided Riemann-Liouville fractional derivative of order a is defined as 

^^£>f JV(a;, t) = [±f{I^-<^N{x, t)), ^{f,) > 0, n = [SR(/x)] + 1 (A4) 

where [x] represents the greatest integer in the real number x. Caputo derivative (Caputo, 1969) is defined 
in the form 

^D^f{x,t) = / f^""Kx,t) ^ _ 1 < i)fj(^) ^ ^ g ^ (^5) 

° * ■' ^ ' ^ r(m - a) Jo [t - t)«+i-'" V / > \ j 

Qtm ' ^ ' 

where ^^f{x, t) is the m-th partial derivative of f{x, t) with respect to t. When there is no confusion, then 
the Caputo operator ^Df will be simply denoted by oDf . 

A generalization of the Riemann-Liouville fractional derivative operator (A4) as well as Caputo fractional 
derivative operator (A5) is given by Hilfer (2000) by introducing a left-sided fractional derivative operator 
of two parameters of order < < 1 and type < < 1 in the form 



oD^;''iV(x,f) 



(A7) 
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For u = 0, (A7) reduces to the classical Ricmann-Liouvillc fractional derivative operator. On the other hand 
for = 1 it yields the Caputo fractional derivative operator defined by (A5). The Laplace transform formula 
for this operator is given by Hilfer (2000). 

L[oD^-:^'N{x, t); s] = s''iV(a;, s) - s'^^^-^h^^-"^^^-"^ N{x, 0+), < < 1 {AS) 

where the initial value term /q^ '^^^^ '^''iV(a;, 0+) involves the Riemann-Liouville fractional integral operator 
of order (1 — !/)(!— /x) evaluated in the limit as f — > 0+, it being understood that the integral 

N{x,s) = L{N{x,ty,s} = e-''^N{x,t)dt, {A9) 

Jo 

where 5i(s) > 0, exists. 

Note Al. The derivative defined by (A7) also occurs in recent papers by Hilfer (2003, 2009), Srivastava 

et al.(2009), Tomovski et al. (2011) and Saxcna ct al. (2010). 

Following Feller (1952, 1971), it is conventional to define the Riesz-Feller space fractional derivative of 
order a and skewness in terms of its Fourier transform as 

FUD'^fix); k} = -i,i{k)r{k), (AlO) 

where 

^l}l{k) = |A;|"exp[i(si5n k)—], < a < 2. |6'| < min{a,2 - a). {All) 
When ^ = we have a symmetric operator with respect to x, that can be interpreted as 

^D'S = -[-^\-- {A12) 

This can be formally deduced by writing — (fc)" = — (fc^)'^. For 6 ^ we also have 

FUD'^^f{x);k] = -\krr{k). (A13) 

For < a < 2 and \6\ < min{a, 2 — a} the Riesz-Feller derivative can be shown to possess the following 
integral representation in x domain: 



For = 0, the Riesz-Feller fractional derivative becomes the Ricsz fractional derivative of order a for 
1 < a < 2 defined by analytic continuation in the whole range 0<a<2,a^l, see Gorenflo and Mainardi 
(1999), as 

= -A[J+" - II"] {AU) 

where 

A = - ^— -r; -ri" = -^/|""- (^15) 

2cos(a7r/2) =^ da;^ =^ ^ ^ 

The Weyl fractional integral operators arc defined in the monograph by Samko et al. (1990) as 

{4n){x) = fjx - 0^-'N{Od^, /3 > 0; (^16) 

{I^-N){x) = - x)f'-'N{Od^, /S > 0. {A17) 

Note A2. We note that is a pseudo differential operator. In particular, we have 

^Dl = ^, but + A. (A18) 
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